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ABSTRACT 

We have previously identified 10 M31 black hole candidates (BHCs) in M31, from their X-ray 
properties alone. They exhibit "hard state" emission spectra that are seen at luminosities <10% 
Eddington in X-ray binaries (XBs) containing a neutron star (NS) or black hole (BH), at luminosities 
that significantly exceed the NS threshold. Nine of these are associated with globular clusters (GCs); 
hence, these are most likely low mass X-ray binaries (LMXBs); eight are included in this survey. We 
have recently discovered that analysis of the long term 0.5-4.5 keV variability of XBs via structure 
functions allows us to separate XBs from AGN, even though the emission spectra are often similar; 
this has enabled us to search for BHCs outside of GCs. We have identified 26 new BHCs (12 strong, 
14 plausible) within 20' of the M31 nucleus (M31*), using 152 Chandra observations spaced over 
~13 years; some of our classifications were enhanced with XMM-Newton observations. Of these, 7 
appear within 100" of M31*; this supports the theory suggesting that this region experiences enhanced 
XB production via dynamical processes similar to those seen in GCs. We have found a parameter 
space where our black hole candidates are separated from Galactic neutron star binaries: we show 
that modelling a simulated hard state spectrum with a disk blackbody + blackbody model yields 
parameters that lie outside the space occupied by neutron star binaries that are modeled this way. 
The probability that our BHCs all lie within the NS parameter space is ~ 3 x 10~ 29 . 

Subject headings: x-rays: general — x-rays: binaries — black hole physics 



^ _ 1. INTRODUCTION 

We know of ^20 X-ray binaries with dynamically con- 
firmed black hole (BH) accretors; these include 15 low 
mass X-ray binaries (LMXBs) and 3 high mass X-ray 
binaries (HMXB s) in the Milky way and Mag ellanic 
^J-" ; Clouds (see e.g. iRemillard k McClintockl 120061 and 
f^S . references w ithin), as well as BH + Wo lf-Rayet bina- 
Cn ■ ries in IC10 (iSilverman k Filippenkol l2008l) and NGC300 
t— I \ (|Crowther et al-lfeOlOl ) . These systems were identified us- 
ing X-ray and optical observations; the mass function is 
. £h ' calculated from periodic radial velocity shifts in emission 
S^ lines from the optical counterpart. All of the BH LMXBs 
identified this way are necessarily transient, because the 
C3 ■ optical spectra of bri ght LMXBs are domina ted by th e 
accretion disc (see e.g. van Paradijs k McClintock 1994). 
We have established a method for identifying BHCs 
from their X-ray properties alone. This makes use of 
the "low/hard" emission state seen in BH and neutron 
star (NS) XBs (|van der Klisl 119941 ) . that is only seen 
at 0. 01-1000 keV luminosit i es <10% Eddingt on in NS 
XBs (jGladstone et al.l[2007l ): rTang et al.l (f201ll ) recently 
found that the low/hard state is limited to luminosi- 
ties <10% Eddington in BH XBs also. BH XBs may 
exhibit low/hard emission states at considerably higher 
luminosities than NS XBs, due to the higher accretor 
mass. However, it is necessary to differentiate between 
our BHCs and distant active galactic nuclei (AGN), since 
AGN and XB emission spectra are often similar. 

We have identified 10 BHCs from their high lumi- 
nosity low states to date. Of these, 9 are associated 



with M31 GCs (IBarnard et all [2001 IBarnard k Kolbl 
I2009t IBarnard et al.l 1201 lal l2012bl) . and are therefore 
likely LMXBs. Only one of these is transient; how- 
ever, persistently bright GC BH XBs are consistent with 
tidal capture of a main sequence donor (Kalogcr a et al.l 
120041 although the donor may be disrupted in the pro- 
cess), or w ith an ultra-compac t system with a degener- 
ate donor (jlvanova et al.|[2~010D . We identified our first 
BHC outside of a GC from its X-ray spectra, long term 
(~12 year) behaviour, and a serendipitous HST observa- 
tion (jBarnard et alJ[2011b[ ) ; the faint optical counterpart 
(Mb > —0.4) suggests a low mass donor for this system 
also. We cover 8 of the 9 GC BHCs in our survey: BHCs 
1, 2, 20, 25, 28, 3 1 , 32, a nd 34. The field BHC described 
in IBarnard et al.l (|2011bf ) is BHC3. 

1.1. Our GC BHCs 

We obtained accurate positions for our X-ray sources 
by registering 27 X-ray sources associated with M31 
GCs to the M31 F ield 5 B band image provided by 
iMassev et all (J2006I) . The r.m.s. offsets w ere 0.11" in 
R.A. and 0.09" in Dec (jBarnard et al.ll2012al ): this would 
be extremely unlikely unless the 27 X-ray sources used 
for calibration were indeed associated with the GCs. 

The X-ray spectra of BH binaries are usually described 
with two components: a thermal component (often mod- 
elled as a multi-temperature disk blackbody), and a 
power law component to represent unsaturated, inverse- 
Com pton scattering of cool p hotons on hot electrons 
(|Remillard k McClintockl l2006f ). The hard state is clas- 
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sified by a power law component with photon index (T) 
= 1.4-2.1, and a thermal co mponent that contributes 
< 20 % of the 2-20 keV flux (jRemillard fc McClintockl 
[200fl . 

Our very first GC BHC (XB045, XBo 45 in 
iBarnard et al.l [20081 ) was associated with the M31 GC 
B045, name d following the Rev i sed Bologna Cata logue 
v.3.4 (RBC, IGalletiet al.ll200l 120061. [20071 [2009h . Its 
^17,000 count XMM-Newton/pn spectrum was well de- 
scribed by an absorbed power law with line-of-sight ab- 
sorption (N H ) = 1.41±0.11xl0 21 atom cm -2 , and T = 
1.45±0.04; x 2 /dof = 517/487 (good fit probability, gfp, 
= 0.17). Adding a blackbody component improved the 
fit somewhat, but the thermal contribution to the flux 
was too small to be constrained; we therefore considered 
XBo 45 to be in its hard state. Since its unabsorbed 
0.3-10 keV luminosity was 2.5±0.2xl0 38 erg s" 1 , 140% 
Eddington for a 1.4 M Q NS, we considered it a BHC. 

The 8 remaining GCs BHCs are included in this sur- 
vey, unlike XBo 45, so we don't describe each spectrum in 
detail here. However, we note th at fitting a two compo- 
nent model to XB144 (XBo 144 in lBarnard fc Kolbll2009ft 
yielded kT = 0.0082±0.0016 keV, indicating a complete 
lack of thermal component in the spectrum. We also 
note that even though XB082 was best described by 
T = 1.20±0.09, the x 2 /dof for that fit was ~0.9, and 
we were able to obtai n fits where T=1.4 and x 2 /dof <1 
([Barnard et al.ll2011al ). 

The unabsorbed 0.3-10 keV luminosities of our GC 
BHCs exhibiting hard state spectra r a nge over ~5- 
45xl0 37 erg s" 1 ([Barnard fc Kolfj 12001 IBarnard et ahl 
I2011al I2012bl ). Comparison with the 0.5-10 keV AGN 
flux distribution obtained by iGcorgaka kis et all ([20081 ) 
yields a 2.4xl0~ 36 probability that our GC BHCs are 
coincident AGN. The probability that our brightest GC 
BHC, XB135, is a coincident AGN is 1.2xl0~ 6 . We will 
present evidence in a separate paper that XB135 may 
contain the most massive stellar mass BH known to date; 
it may have been formed by direct collapse of a high mass, 
metal poor star (R. Barnard et al., 2013, in prep). 

We found that the GCs hosting these very bright X- 
ray sources were significantly more massive and/or metal 
rich than the other GCs in M31, agreeing with previous 
work. However, two GCs had particularly low metalici- 
ties; one of these is B135, consistent with the direct col- 
lapse f ormation scenario for X B135 (see IBarnard et ahl 
I2012U and references within). [Bclczvns ki. et alT (2010) 
have shown that while BH masses are limited to ^15 M® 
for Solar metalicties, they could theoretically reach ~80 
Mq for metalicities ^0.01 Solar. 

1.2. Could neutron star binaries mimic the BH hard 

state? 

High accretion rate NS XBs exhibit multi-component 
emission that may appear to be hard state spectra in 
extragalactic X-ray sources, where the spectra have rel- 
atively few photons. To compare the emission spectra 
of our BHCs with NS XRB, we use the double thermal 
(disk blackbody + blackbody) model used by Lin et al. 
(2007, 2009, 2012) to describe hundreds of RXTE obser- 
vations of NS XBs in all their varied emission states. We 
note that the physical interpretation of their model is 
contradicted for persistent XBs by a substantial body of 



work; however, Lin et al. (2007, 2009, 2012) have sam- 
pled spectra from the full range of NS LMXB emission 
states in a consistent manner, allowing us to compare our 
BHCs with NS XBs in a single parameter space. 

For a long time, NS XBs were divided into two types, 
Z-sources (high lumonisity) and atoll sources (low lumi- 
nosity), based on their luminosities and color-color dia- 
grams (CDs); the CDs of Z-sources had three branches 
(horizontal, normal, and flaring), and evolved along these 
branches without ever jumping from one branch to the 
other; the CDs of the atol l sourc es were more fragmented 
(jHasinger fc van der Klisl 119 89) . Furthermore, the Z- 
sources were split into those like Cygnus X-2, an d those 
like Scorpius X-l ([Hasinger fc van der Klisl[l989|) . It was 
believed that the differences were due to more than just 
the accretion rate, because the Z-sources varied by a fac- 
tor of a few when tracing their Z-shaped CDs, while 
atoll source intensities varied by 1-2 orders of magni- 
tude (|Muno. Remillard. fc Chakrabartvl2002fl . However, 
we now know of two transient NS systems that exhib- 
ited both types of Z-source behavio r before evolving to 
atoll source behavior du ring decay (Ho man et al.l [20071 : 
Chakraborty. Bhattac haryya. fc M ukhcriec 20111 ). 

iLin et all (|2007l ) examined the spectral evolution of 
two Galactic X-ray transients, Aql X-l and 4U 1608—52, 
over many RXTE observations covering >20 outbursts. 
They devised a new double-thermal model (disk black- 
body + blackbody) to describe a NS transient soft 
state that is analagous to the BH so ft state de- 
scribed by |RemiIIard!&!McCh^oc3 $20M)- They have 
since applied their model to RXTE observations of 
XTE J1701-462, one of the transients that exhibits Cyg- 
like and Sco-like Z- source behavior as well as atoll behav- 
ior (ILin et al.ll2 009). and also to the Sco-like Z-source GC 
17+2 (Li n et ahl 12012( 1. They have applied their model 
to hundreds of RXTE spectra from Galactic NS binaries 
including the full gamut of NS spectral behavior. 

They found that their disk blackbody + blackbody was 
unsuccessful in two situations. Firstly, they found the 
hard st ate spectra to be pow er law dominated, as ex- 
pected ([Lin et al.ll2007L I2009T ). We therefore expect our 
BHCs to inhabit a separate parameter space to the NS 
XBs because fitting the double thermal model to hard 
state spectra will yield unphysical results. Secondly, they 
found that Z-sources required a three-component spec- 
trum (disk blackb ody + blac kbody + po wer law) on the 
h orizontal branch Jlin et al.ll2009Ll2012H . 

ILin et ahl (J20101 ) also performed broadband analysis of 
Suzaku and BeppoSAX observations of the persistently 
bright Galactic NS XB 4U 1705—44, with energy ranges 
0.1-600 keV and 0.1-300 keV respectively. When using 
a disk blackbody + blackbody model to the soft state 
spectra, they found temperatures that were similar to 
those they observ ed in the RXTE observations ([Lin et al.1 
120071 120091 120121) . Fitting an additional Comptonization 
component, via eith er a power law or SIMPL convolu- 
tion model ffollowing lSteiner et a l. 2009) resulted only in 
small changes in the thermal components; this is because 
the Comptonized component only contributed 5,10% of 
the t otal flux in the soft state (see Fig. 6 in ILin et all 
l2010f h Hence the parametric differences between our 
BHCs and the NS XBs should not be due to differences 
in energy bands used in the observations, or the lack of 
a third (power law) component. 
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Evidence against the double thermal model indicates 
an extended corona that contributes a substantial por- 
tion of the X-ray flux. This evidence includes the ingress 
times of periodic absorption dips in the X-ray lightcurves 
of the high in c lination XBs (the dipping sources, s ee e.g. 
IChurchl l200lt IChurch fc Balucihska-Churchl |20M and 
references within), and also broadened emi ssion lines in 
a Cha ndra grating spectroscopy of Cyg X-2 ()Schulz et al.l 
2009). Furthermore, compact corona models where the 
inner disk temperature is tied to the seed photon energy 
for Comptonization are rejected for ULXs in NGC253, 
and the confir med BH+Wolf-Rayet binary IC10 X-l 
(Barnard 2010) , as well as B HC3 in the steep power law 
state (jBarnard et al.|[20lTb] ). 



1.3. Structure function analysis of GC XBs 

In IBarnard et al.l (|2012b[ ). we applied structure func- 
tion analysis to XB s for the first time, following 
iVagnetti et al.l (120111 ). who created an ensemble struc- 
ture function for AGN. They used a structure function 
(SF) to estimate the mean intensity deviation for data 
separated by time t: 



5F(r) s y|(|log/ x (t + r)-log/ x (i)|) 2 -a2, (1) 

where o~ n is the photon noise and fx is the X-ray flux. 
They grouped the SF into logarithmic bins with width 
0.5; each bin in the range log(r) = 0.0-3.0 contained 
more than 100 measurements. 

Our sample consisted of 37 X-ray sources associ- 
ated with objects in the RBC; these were classified by 
ICaldwell et all ((20091) as 30 confirmed GCs, 4 GC can- 
didates, 1 star, and 2 AGN. The SFs of GC XBs with 
0.3-10 keV luminosities ~2-50xl0 36 erg s _1 tended to 
show significantly more variability than AGN over a wide 
range of time-scales. The SFs of brighter XBs gener- 
ally showed comparable or less variability than AGN, 
despite their high signal to noise; however, their high 
fluxes made them unlikely AGN. Hence, SFs provide an 
effective mechanism for distinguishing between XBs and 
AGN (|Barnardet al.ll2012bl) . 

1.4. Searching for field BH XBs 

In this work, we combine these techniques to search 
for BHCs in the whole region covered by 152 Chandra 
observations from our monitoring programme; the roll 
angle is unrestricted, resulting an approximately circular 
field of view with radius ~ 20'. We surveyed 530 X- 
ray sources in this region. We also obtain spectra from 
archival XMM-Newton observations, to strengthen the 
cases for several BHCs. 

The 100" region surrounding M3 1 * is p articularly in- 
teresting, because IVoss fc Gilfanovl (|2007t ) found an ex- 
cess of X-ray binaries over the radial distribution ex- 
pected from K band light (tracing stellar mass); this ex- 
cess h ad the distribution expected of dynamically formed 
XBs (Fabian e t alj[l975l) . Dynamical XB formation re- 
quires stellar densities rarely seen outside GCs, and 
IVoss fe Gilfanovl (J2007I) proposed that the M31 bulge is 
sufficiently large and dense to form a significant number 
of XBs dynamically. However, since the stellar velocities 
in the M31 bulg e are ~ 5-10 times higher than in GCs, 
Voss & Gilfanov (2007) expect only short period binaries 



to survive, with most of these containing BH accretors. 
Evidence for BHCs in this region would therefore provide 
strong support for their theory. 

In the next section, we discuss the observations and 
data analysis; this is followed by our results, then by our 
discussion. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Chandra analysis 

The central region of M31 has been observed with 
Chandra on a ~ monthly basis for the last ~13 years 
in order to monitor transients; we exclude periods when 
M31 cannot be observed due to satellite constraints (~ 
March- April) . We have analyzed 98 ACIS observations 
and 54 HRC observations. We determined the position 
of each source from a merged 0.3-7.0 keV ACIS image, 
using the IRAF tool imcentroid. The X-ray positions 
were registered to the Field 5 B band imag e of M31 from 
the L ocal group Galaxy Survey (LGS, see iMassev et al.l 
120061) , using 27 GC X-ray sources a n d follow ing the pro- 
cedure described in IBarnard et al.l (|2012al ). The r.m.s 
uncertainties in registration were 0.11" in R.A. and 0.09" 
in Dec. 

We obtained 0.3-7.0 keV lightcurves and spectra from 
circular source and background regions for each source. 
The background region was the same size as the source 
region, and at a similar off-axis angle. The extraction 
radius varied between sources, because larger off-axis 
angles resulted in larger point spread functions. We 
used the CIAO v4.5 software suite, with corresponding 
CALDB to reduce the data, and XSPEC vl2.7 to analyze 
the spectra. 

For ACIS observations, we obtained corresponding re- 
sponse matrices and ancillary response files. We initially 
estimated the conversion from flux to luminosity by as- 
suming a power law emission spectrum with photon in- 
dex 1.7, with A^h = 7xl0 20 atom cm -1 , then determin- 
ing the unabsorbed 0.3-10 keV luminosity equivalent to 1 
count s _1 at the location of the source. After correcting 
for the exposure, vignetting and background, multiply- 
ing the source intensity by this conversion factor gave the 
source luminosity. Source spectra with >200 net counts 
were freely fitted. When good spectral fits were found 
for a particular observation of a source, the parameters 
of these fits replaced the default parameters when es- 
timating the source luminosity for observations of the 
source with <200 photons. 

For HRC observations, we we included only PI chan- 
nels 48-293, thereby reducing the instrumental back- 
ground. We used the WebPIMMS tool to find the un- 
absorbed luminosity equivalent to 1 count s , assum- 
ing the same emission model as for the ACIS observa- 
tions with <200 photons. We created a 1 keV exposure 
map for each observation, and compared the exposure 
within the source region with that of an identical on-axis 
region, in order to estimate the estimate the necessary 
exposure correction. We multiplied the background sub- 
tracted source intensity by the correction factor to get 
the 0.3-10 keV luminosity. 

We created long term 0.3-10 keV lightcurves for each 
source, using the luminosities obtained from each obser- 
vation as described above; all luminosities assume a dis- 
tance of 780 kpc (|Stanek fc Garnavi ch 1998). We only 
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Table 1 



For each BHC we provide its position, and the number of ACIS and HRC observations (Oa and Oh respectively). We then give the best 

fit constant 0.3—10 keV luminosity, normalized to 10 37 erg s _1 , with the corresponding x 2 /dof. We finally show the number of 

observation pairs used to create the structure function. Uncertainties are quoted at the ltr level. 



BHC 


Position 


Oa 


O h 


Lbp/IO 3 ' 


xVdoi 


7V SP 


1 


00:42:15.786+41:01:14.24 


18 


48 


21.3+0.4 


191/65 


2145 


2 


00:42:18.648+41:14:01.85 


83 


51 


5.54+0.06 


1479/136 


9316 


3 


00:42:22.954+41:15:35.23 


91 


51 


7.57+0.06 


23850/144 


10440 


•1 


00:42:26.070+41:19:15.03 


93 


51 


2.53+0.04 


3576/146 


10731 


5 


00:42:28.193+41:10:00.30 


85 


51 


3.85+0.05 


486/138 


9591 


6 


00:42:28.285+41:12:22.95 


89 


51 


4.73+0.05 


6113/142 


10153 


7 


00:42:31.147+41:16:21.67 


93 


51 


9.32+0.08 


1123/146 


10585 


8 


00:42:34.669+40:57:14.20 


13 


12 


0.64+0.08 


696/24 


300 


9 


00:42:39.585+41:16:14.30 


70 


21 


0.063+0.006 


3480/90 


4095 


10 


00:42:40.654+41:13:27.32 


93 


53 


1.40+0.03 


4724/145 


10585 


11 


00:42:42.177+41:16:08.23 


74 


31 


0.057+0.006 


2331/107 


5778 


12 


00:42:44.831+41:11:37.89 


92 


51 


3.6+0.04 


3220/145 


10585 


13 


00:42:45.122+41:16:21.68 


96 


51 


2.3+0.03 


3291/149 


11175 


11 


00:42:45.946+41:10:36.53 


51 


27 


0.035+0.007 


1759/80 


3240 


15 


00:42:46.969+41:16:15.58 


95 


51 


3.1+0.04 


3202/148 


11026 


1(5 


00:42:47.176+41:16:28.41 


92 


51 


0.294+0.014 


40403/145 


10585 


IV 


00:42:47.893+41:15:32.87 


95 


51 


2.6+0.03 


3278/148 


11026 


18 


00:42:48.529+41:15:21.12 


95 


51 


11.6+0.09 


2369/148 


11026 


19 


00:42:48.546+41:25:22.12 


60 


51 


1.6+0.05 


8701/110 


6105 


20 


00:42:52.030+41:31:07.87 


11 


54 


49.4+0.4 


343/64 


2080 


21 


00:42:52.534+41:18:54.46 


95 


51 


11.79+0.09 


1271/148 


11026 


22 


00:42:54.935+41:16:03.19 


95 


51 


4.8+0.05 


18769/148 


11026 


23 


00:42:57.900+41:11:04.65 


93 


51 


5.3+0.05 


826/146 


10731 


21 


00:42:59.675+41:19:19.35 


93 


51 


5.4+0.05 


887/146 


10731 


25 


00:42:59.872+41:16:05.72 


95 


51 


5.6+0.06 


1074/148 


11026 


26 


00:43:02.937+41:15:22.53 


95 


51 


3.3+0.04 


4838/148 


11026 


27 


00:43:03.220+41:15:27.69 


95 


51 


1.5+0.03 


11960/148 


11026 


28 


00:43:03.876+41:18:04.91 


94 


51 


4.1+0.05 


2920/147 


10878 


29 


00:43:05.667+41:17:02.43 


77 


29 


0.053+0.005 


10972/105 


4851 


30 


00:43:09.866+41:19:00.76 


85 


18 


0.127+0.010 


2675/132 


8646 


31 


00:43:10.622+41:14:51.30 


90 


54 


13.4+0.11 


915/143 


10296 


32 


00:43:17.593+41:27:44.87 


17 





0.11+0.02 


409/16 


136 


33 


00:43:34.334+41:13:23.08 


45 


54 


3.8+0.06 


1380/98 


4851 


31 


00:43:37.320+41:14:43.09 


34 


54 


6.7+0.09 


209/87 


3828 


35 


00:43:44.593+41:24:10.25 


32 


39 


0.5+0.03 


2296/71 


2346 



included observations with net source counts > af- 
ter background subtraction. We fitted each long term 
lightcurve with a line of constant intensity, in order to 
ascertain the source variability. 

We note that the ratio of HRC to ACIS luminosity de- 
pends strongly on the spectral model, and should be 1.0 
if the model is correct. Differences between actual and 
assumed emission spectra during HRC and faint ACIS 
observations may lead to systematic offsets between lu- 
minosities from adjacent HRC and ACIS observations. 

We derived SFs from the 0.5-4.5 keV fluxes of each ob- 
servation of every target by assuming a power law spec- 
trum with the same photon index as for the HRC and 
faint ACIS observations; an M31 X-ray source with a 
0.3-10 keV unabsorbed luminosity of l.OxlO 37 erg s _1 , 
A H = 7xl0 20 atom cm" 2 and T = 1.7 has a 0.5-4.5 keV 
fl ux of 0.80 xlO" 1 3 erg cm" 2 s _1 . 

iVagnetti et al.1 (|2011h calculated the noise component 
from 



(2) 



our case. 



ai = 2 ( ( 51og/x) z )^2(loge) 




n0.5 



assuming that Sfx/fx = (l/A pllo t) '° , and N phot is the 
number of photons. Our lightcurves are background sub- 
tracted, and ARF-corrected; furthermore, uncertainties 
in the luminosities of bright sources include uncertainties 
in the spectral parameters. As a result, our uncertainties 
are not simply due to photon counting noise. Hence in 



(log e) 



afx (t + t) 



fx (t + t) 



<rfx (*) 



fx(t) 



(3) 



where afx is the uncertainty in X-ray flux and often sig- 
nificantly larger than N®^ ot . In the simplest case where 
only 1 pair has a particular separation, then 



(loge)' 



ill 
h 



oh 
h 



(4) 



where /i and fi are the fluxes of the two observations. 

2.2. XMM-Newton analysis 

In addition to our Chandra observations, we also ob- 
tained spectra for some of our BHCs from archival XMM- 
Newton observations. Their higher sensitivity and longer 
exposure times often yielded superior spectra for bright 
sources. We used only the pn instrument, because the 
MOS detectors are more prone to pile up. We used SAS 
vlO for the data reduction. 

We filtered out background flares by cre- 
ating a lightcurve with the selection "(PAT- 
TERN==0)&&;(FLAG==0)&&(PI in [10000:12000])" 
and 100 s binning, then excluding intervals with 
>0.4 count s^ 1 . Source and background spectra were 
obtained from circular regions with 20-40" radius, 
along with corresponding response files; the filter 
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Unfolded Spectrum 
kT (DBB) = 0.67 keV, kT(BB) = 1.76 keV F -0.S 




1 2 

Energy (keV) 

Figure 1. Faked spectrum created using an absorbed power law 
model that represents the hard state: JVjj = 7xl0 20 atom cm -2 , 
r = 1.7. The best fit double thermal model gives a 0.67 keV disk 
blackbody and a 1.76 keV blackbody; the blackbody component 
contributes ~80% of the 2-10 keV flux. 

"(PATTERN<=4)&&(FLAG==0)" was applied along 
with the good time filter and extraction region. 

2.3. Identifying black hole candidates 

Our BHC identification process involved three steps. 
First, we had to establish that the X-ray source was an 
XB rather than an AGN. Then we checked the ACIS 
observations of each XB, to see if it exhibited a high lu- 
minosity hard state (r < 2.1 at a 0.3-10 keV unabsorbed 
luminosity >3xl0 37 erg s _1 , i.e. 10% Eddington for a 2 
Mq NS). Finally we found the highest quality observa- 
tion that exhibited a high luminosity hard state, fitted 
it with the double-thermal model, and compared the re- 
sults with the known range of parameters exhibited by 
the NS systems studied by Lin et al. (2007, 2009, 2012). 

Likely XBs were identified from their associations with 
GCs, variable SFs or relati vely high fluxes. Th ere are 
-170 GCs within our field (IPeacock et. al.ll2011l): using 
the .5-10 keV AGN flux distribution (Gcorgakakis et al. 
2008), we expect to find 0.0005 AGN with flux 



>1.4xl0~ 13 erg cm 2 s ' 



(10 37 erg s" 1 ) within 1" of 



any GC. Hence any bright X-ray source associated with 
a GC is a likely XB. We also identify XBs using SFs 
that are strikingly more variable than the ensemble AGN 
SF. Furthermore, only — 1 AGN per square degree is ex- 
pected with a flux equivalent to ~ 5 x 10 37 erg s _1 ; hence 
brighter X-ray sources are likely XBs too. 

Lin et al. (2007, 2009) found that they were unable to 
successfully fit hard state spectra with their double ther- 
mal model. This is because the disk blackbody compo- 
nent must account for the low-energy flux, resulting in an 
unphysically low temperature. To demonstrate this, we 
created a fake hard state spectrum using the XSPEC com- 
mand fakeit; we used a power law with T = 1.7, with 
absorption equivalent to 7xl0 20 atom cm~ 2 . We present 
the best fit double thermal model for this spectrum in 
Fig. [TJ which consists of a 0.67 keV disk blackbody and 
a 1.76 keV blackbody; the blackbody component con- 
tributes -80% of the 2-10 keV flux. By contrast, the 
soft states of NS XBs studied by Lin et al. (2007, 2009, 
2012) yielded disk blackbody temperatures >1 keV, with 
the blackbody component contributing less than 50% of 
the 2-10 keV flux. 



Since Lin et al. (2007, 2009, 2012) failed to fit hard 
state spectra with the double thermal model, we ex- 
pect disk blackbody + blackbody fits to our BHCs to 
inhabit a different parameter space from the NS binaries 
if they really are in the hard state. If a BHC exhibits a 
blackbody temperature <1.5 keV, a disk blackbody tem- 
perature <1.0 keV (<1.2 keV for disk blackbody 2-10 
keV luminosity >2xl0 37 erg s _1 ), and/or a blackbody 
contribution to the total 2-10 keV flux >45% at a 3cr 
level, then we consider it a strong BHC; otherwise, we 
label it a plausible BHC. These criteria are drawn from 
the hundreds of RXTE spectra of Aql X-l, 4U 1608-52 
XTEJ1701-462, and GX 17+2 analyzed by iLin et all 
(|2!Ml2tM 12012 ) . 

3. RESULTS 

TableQ]gives the location of each of our BHCs, followed 
by the number of observations with ACIS and HRC. Next 
we provide the best fit line of constant 0.3-10 keV lu- 
minosity over —13 years, with the corresponding \ 2 and 
number of degrees of freedom (dof). Finally, we show the 
number of observation pairs used to create the structure 
function. Usually this is -/V bs (-Nobs — -0 /2j where A b s 
is the total number of observations; however, BHCs 7, 
29, 30, and 35 included observations with zero luminos- 
ity that could not be included in the SF. All uncertainties 
in this work are quoted at the la level unless specified 
otherwise. 

For nearly all BHCs, closely-spaced ACIS and HRC 
observations give consistent luminosities, although some 
small systematic offsets occur due to differences between 
assumed and actual spectra when fitting is not possible. 
However, BHC32 appears to have substantial signal in 
the HRC observations even when this transient source 
is in quiescence; the HRC is considerably more sensitive 
than ACIS to low energy photons, hence the HRC obse- 
vations of BHC32 must be contaminated by soft X-rays. 
BHC32 is located near 2 variable stars (16" and 17" dis- 
tant); since the extraction radius for this souce is 15" 
due to the high off-axis angle, our HRC observations are 
likely contaminated by one or both foreground stars. As 
a result, we excluded all HRC observations of BHC32 
from our analysis. 

3.1. BHC locations 

Figure [2] shows a detail of the B band image of M31 
Field 5 from the LGS (Massey et al., 2006), superposed 
with the positions of our BHCs; north is up, east is left. 
The circle encloses the region within 100" (—380 pc) of 
M31*; this region is enlarged in the top right portion of 
the figure. BHCs 1 and 8 lie below the southern edge of 
Field 5. 

3.1.1. The central 100" region 

We find 7 BHCs within 100" of M31*, 20% of our 
total sample: BHCs 9, 11, 13, 15, 16 , 17, and 18. 
This res ult appears to be cons istent with the predic- 
tions of IVoss fc Gilfanovi (|2007j) . BHC9 is a recurrent 
transient with a peak 0.3-10 keV luminosity of — 8xl0 37 
erg s" 1 . BHC11 was bright for —1-2 years, then turned 
off. BHC18 appears to be persistently bright (— 10 38 
erg s _1 ), and may vary rapidly if the HRC luminosities 
are correct; however, the emission from BHC18 is vari- 
able (r —1.4-2.4 for bright ACIS observations), meaning 



Barnard et al. 




Figure 2. Detail of the B band image of M31 Field 5 from the LGS {Massev et al. 2006), superposed with the X-ray positions of our 
BHCs. The blue circle has a 100" radius (~380 pc), and denotes the region where Voss et al. (2007) predict enhancement for dynamically 
formed XBs, many of which are expected to contain black holes. North is up, and East is left. A magnified view of the central 100" region 
is presented at the top right of the image. BHCs 1 and 8 are located outside the Field 5 image. 

that the HRC luminosities may be systematically offset 
due to differences between assumed and true spectra. 

BHCs 13, 15, 16, and 17 are extremely variable, 
with BHCs 13, 15, and 17 luminosities ranging over 
~l-6xl0 37 erg s _1 , and BHC16 exhibiting luminosities 
from ^5xl0 36 to ~2xl0 38 erg s _1 . We observed sim- 
ilar behavior in the X -ray source associated with B158 
(|Barnardet~aLll2012bD, a high inclinati on LMXB with 
a ~10000 s period dTrudolyubov etHI I2002D. a nd an 
asymmetric, precessing disc (|Barnard et al.l 120061 ). Such 
behavior is seen in low mass ratio (short period) sys- 
tems where the outer disc reaches the 3:1 resonance with 
the donor star, causing additional tidal torques that 
lead to elongation and precession of t he disc ( see e.g. 



p 
Osapl989; W hitehurst fc KindllQflTI : IQgilvie fc Dubusl 
B158 exhibited rapid variation over 4-40xl0 37 
erg s , which we speculated might be due to vary- 
ing accretion rates over the disc precession cycle; we 
suggested that other XBs exhibiting such behaviour 
may also be in precessing (i.e. short period) systems 



(|Barnard et al.ll2~012bf ) . Therefore we propose that BHCs 
13, 15, 16, and 17 are shor t period BH XBs, as predicted 
bv lVoss fc Gilfanovi (pOOl . 

3.1.2. Further GC associations 

Two of our new BHCs are associated with confirmed 
old GCs; BHC4 is associated with B096 (following 
the RBC nami ng con vention), and BHC24 with B143. 
iCaldwell et~aTI (J2011I ) rank B096 as the 36th most mas- 
sive and 34th most metal rich out of 379 GCs ([Fe/H] = 
—0.28); B143 is ranked 52nd most massive and 17th most 
metal rich, with [Fe/H] = —0.07. We have found that the 
GCs hosting our previous BHCs were particularly mas- 
sive, or metal rich, or both (Barnard et al., 2011a, 2012); 
hence it is encouraging that the new GC BHCs are asso- 
ciated with GCs that are amongst the most massive and 
most metal rich of the M31 GCs studied by Caldwell et 
al. ( 2009, 2011). 

In Barna rd" et al.l (|2012bT > we classified the semi-regular 
transient associated with HII rich GC candidate SK059A 
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Figure 3. The top half of each panel shows the 0.3-10 keV luminosity of in each Chandra observation made over ~13 years; the red circles 
and grey crosses represent ACIS and HRC observations, respectively. The bottom panel shows T for the best fits to freely- fitted spectra 
with a power law emission model; the mean T is indicated by a line. Uncertainties are at the la level for luminosities, and 90% confidence 
level for i\ An extended, color version of this figure is available in the electronic edition. 



as a possible HMXB associated with the HII, because 
the intervals between its outbursts were consistent with 
n cycles of 120 days; hence the accretor could make an 
eccentric, ^120 day orbit around the high mass donor, 
with mass transfer restricted to times near periastron. 
This system is BHC30. 

3.2. Long term behavior 

We initially identified most of our BHCs from their 
long term (~13 year) light curves and spectral histo- 
ries. They exhibited spectra consistent with the hard 
state (r <2.1) at 0.3-10 keV luminosities >3xl0 37 erg 
s^ 1 in at least one observation. For such spectra, we 



expect the 0.01-1000 keV luminosities to be >3xl0 37 
erg s , our lim i t for n eutron star hard states following 
IGladstone et all (|2007t ). 

We present the ~13 year, corrected, 0.3-10 keV lumi- 
nosity lightcurve for each BHC in Fig. [3l along with its 
spectral history. The top panel in each case shows the 
lightcurve, with ACIS and HRC observations represented 
by circles and crosses respectively. The bottom panel 
shows r for each ACIS observation; sufficiently bright 
observations have freely fitted T, while faint observations 
have r fixed to the mean value; this approach should be 
valid because we expect the faint observations to be in 
the hard state also, with similar a emission spectrum. 
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Table 2 

Details of the highest quality spectrum that indicates a high luminosity hard state. We show the observation number, and net source 

counts. We then show the absorption (-/Vh), spectral index (r), x 2 /dof and 0.3-10 keV luminosity for the best fit power law model. 

Observations starting with zero are XMM-Newton observations; the others are from Chandra. Uncertainties are quoted at the ltr level. 



BHC 


Obs 


Cnts 


JVh/10 2 ^ cm" 2 


r 


X^/doi 


L /10^ Y 


1 


0109270101 


3229 


0.40±0.05 


1.24±0.09 


133/155 [0.89| 27. !-:-!. 1 


2 


303 


704 


0.15±0.10 


1.5±0.2 


33/29 [0.29 


6.6±0.6 


3 


13827 


5597 


0.07 f 


1.50±0.04 


161/182 [0.87] 14.4±0.4 


1 


0690600401 


2370 


0.26±0.04 


1.67±0.14 


153/146 [0.321 6.6±0.3 


5 


13825 


1425 


0.21±0.09 


1.68±0.18 


8/11 [0.72] 


4.2±0.3 


6 


0112570101 


10384 


0.112±0.012 


1.72±0.06 


428/411 [0.26] 7.30±0.17 


7 


0112570101 


14523 


O.llOiO.010 


1.76±0.04 


519/490 [0.181 8.79±0.17 


8 


9524 


352 


0.07 f 


1.29±0.17 


16/15 [0.37 


17±3 


9 


1575 


3097 


0.10±0.03 


1.73±0.09 


106/109 [0.551 5.0±0.2 


10 


13827 


1679 


0.13±0.06 


1.56±0.13 


72/67 [ 


0.32 


3.6±0.2 


11 


303 


704 


0.10±0.09 


1.43±0.18 


58/45 


0.10 


9.7±0.8 


12 


13825 


2620 


0.15±0.04 


2.00±0.10 


90/98 


0.71 


5.2±0.2 


13 


14197 


1593 


0.12±0.05 


1.63±0.14 


72/64 


0.24 


3.2±0.2 


11 


9521 


183 


0.07 f 


1.5±0.4 


6/7 [0.56] 


10±3 


15 


14198 


1860 


0.13±0.05 


1.56±0.12 


75/74 [0.47 


3.9±0.2 


16 


13825 


4003 


0.12±0.03 


1.40±0.07 


161/165 


10.2±0.4 


17 


14197 


1954 


0.14±0.04 


2.05±0.13 


94/76 


3.25±0.18 


18 


13825 


4207 


0.09±0.03 


1.45±0.08 


140/154 


0.78] 9.0±0.4 


19 


0112570101 


5032 


0.41±0.03 


1.89±0.08 


210/226 


0.77] 9.74±0.4 


20 


0402560901 


20025 


0.275±0.009 


1.60±0.02 


771/769 


0.44 45.0±0.4 


21 


13827 


3847 


0.22±0.05 


1.73±0.09 


127/142 


0.82 10.0±0.4 


22 


0112570101 


10145 


0.134±0.014 


1.81±0.05 


490/407 [3 E- 


31 8.2±0.2 


23 


14198 


2508 


0.13±0.04 


1.64±0.10 


111/95 [0.12] 5.1±0.2 


21 


0112570101 


6804 


0.120±0.015 


1.91±0.06 


265/294 [0.89] 5.13±0.15 


25 


0112570101 


5116 


0.108±0.019 


1.50±0.06 


254/236 [0.201 5.3±0.17 


26 


13825 


1425 


0.16±0.07 


1.71±0.17 


54/57 


0.60 


3.3±0.2 


27 


14197 


2099 


0.14±0.04 


1.83±0.12 


69/75 


0.67 


4.1±0.2 


28 


13826 


635 


0.07 f 


1.73±0.13 


23/27 


0.69 


5.3±0.4 


29 


1577 


1375 


0.07 f 


1.33±0.08 


79/58 


0.04 


31±2 


30 


13299 


274 


0.07 f 


1.2±0.2 


10/10 


0.42 


9.9±1.6 


31 


0112570101 


14338 


0.081±0.009 


1.64±0.04 


450/488 [0.891 10.4±0.2 


32 


8184 


521 


0.07 f 


1.43±0.13 


27/23 [0.25 


13.4±1.5 


33 


13825 


1411 


0.14±0.08 


1.81±0.17 


57/60 [0.57 


4.4±0.3 


31 


0112570101 


5181 


0.124±0.019 


1.64±0.07 


246/230 [0.221 7.0±0.2 


35 


14198 


1646 


0.18±0.08 


1.75±0.16 


85/69 [0.10 


5.7±0.4 



BHC13 1 1 175 total pairs, 4560 AC1S pairs 




log(x/day) 

Figure 4. Histogram showing the number of observations with 
separation r for BHC13, which was most frequently observed. The 
histogram with filled circles includes ACIS and HRC observations, 
while the histogram with open circles uses only ACIS observations. 



Figure Q] shows histograms of the number of observa- 
tion pairs with separation r for the most frequently ob- 
served BHC, BHC13. The histogram with filled circles 
includes ACIS and HRC observations, while the open 
circle histogram contains only ACIS observations. The 
observation pairs were logarithmically binned by r, with 
a bin width of 0.2 dex. 

We present the 0.5-4.5 keV structure function for each 
of our BHCs in Fig. [5] with same scheme as for Fig. |4j 
The ensemble SF for typical AGN obtained by Vagnetti 



et al. (2011) is represented by a dashed line in each panel: 
SF(r) = 0.11t 010±001 . We included an ACIS only SF 
for each source in addition to an ACIS + HRC SF be- 
cause faint HRC observations add extremely large noise 
components for some BHCs, suppressing the variation; 
BHCs 8 and 35 only show significant excess variability 
over typical AGN in their ACIS only SFs. 

We find that 28 out of 35 BHCs exhibit SFs with sig- 
nificantly more variation than the Vagnetti et al. (2011) 
ensemble AGN SF (>4cr excess in at least one r chan- 
nel, >3<7 excess in at least 2 channels, or >2cr excess 
in at least 3 channels over the 3<7 limit for the AGN 
SF, assumed to be 0.1 It 013 ). This distinction is im- 
portant, because AGN and XBs often have very similar 
spectra. The X-ray sources with SFs exhibiting similar 
or less variation than the ensemble AGN SF (BHCs 2, 
5, 7, 21, 23, 24, 25, and 34) had mean 0.3-10 keV fluxes 
>5.5xl0 -13 erg cm ~ 2 s _1 . 

We have found 24 X-ray sources in our observations 
that exhibited mean 0.3-10 keV fluxes >5.5xl0 -13 erg 
cm~ 2 s _1 over the ~13 years of Chandra monitoring; 
this includes 8 GC XBs and 1 possible supernova rem- 
nan t. The 0.5-10 keV flux distribution of AGN found 
by (jGeorgakakis et all l2008f ) predicts 0.6 AGN at this 
flux level within our observed region (approximated by 
a circle with radius ~20' or ^4 kpc) . Since the observed 
number of unidentified X-ray sources at these flux lev- 
els is larger than the predicted AGN number by a factor 
~25, we conclude that they are likely XBs. 
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BHCl N H = 5.0+/-0.9 (5/16) r=I.50+/-0.14(5/16) 



BHC2N„ = I.5+/-0.7(5/14)r = 1.43+/-0.11 (6/14) 
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Figure 5. Structure functions (SFs) for each BHC. In each case we provide the absorption (Njj) and photon index (T) used to convert 
from 0.3—10 keV to 0.5-4.5 keV flux, along with the x 2 /dof, obtained by fitting all ACIS observations with >200 counts; Nh is normalized 
to 10 21 H atom cm -2 . The SF with filled circles is drawn from ACIS and HRC observations, while the SF with open circles uses only ACIS 
data. The dashed line represents SF(r) = O.llr °, our approximation of the Vagnetti et al. (2011) ensemble AGN SF. An extended, 
color version of this figure is available online. 



3.3. Spectral analysis 

After the initial identification of candidates, we sought 
the observations that best represented the BHC case for 
each object; these were the highest quality spectra that 
appeared to exhibit high luminosity hard states. Ta- 
ble [2] provides the details of these observations. For 
each BHC we give the observation number, and net 
sources counts. We then give the absorption (Nh), 
photon index (T), x 2 /dof and 0.3-10 keV luminosity 
for the best fit absorbed power law model. Observa- 
tions 303-13827 are Chandra ACIS observations, while 



0109270101, 0112570101, and 0402560901 are XMM- 
Newton observations. 

We note that the spectrum of BHC 22 rejected simple 
absorbed power law models. It is well described by an 
absorbed disk blackbody + power law model with T; n = 



1.97±o.3 9 keV, and T 



2.16^q 19 , with the power law 

-137 



contributing 61±12% of the 8.0±1.1 xlO 37 erg s" 1 0.3- 
10 keV luminosity; it is consistent with either of two 
canonical black hole states: the hard state, or a steep 
power law state similar to GRS 1915+105 (McClintock 
& Remillard, 2006). 
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Figure 6. Comparison of disk blackbody 2-10 keV luminosity vs 
kTi n for the best fit double thermal models to our BHC spectra. 
The shaded region encompasses the parameter space inhabited by 
hundreds of spectral fits for Aql X-l, 4U 1608-52, XTE J1701-462 
and GX 17+2; together, they represent all types of NS LMXB 
behavior (Lin et al. 2007, 2009, 2012). Filled circles represent field 
BHCs, and open circles represent GC BHCs. We note that the disk 
blackbody component contributes only ~l-35% of the total 2—10 
keV flux for our BHCs, but dominates the spectra studied by Lin 
et al. (2007, 2009, 2012). The probability that the BHCs all lie 
within the shaded region is ~3xl0 . 



We applied a disk blackbody + blackbody emission 
model to our BHC spectra. The disk blackbody compo- 
nents were generally better constrained than the black- 
body components, particularly for spectra with relatively 
few net source counts; this is simply because there are 
fewer high energy photons than low energy photons and 
the disk blackbody is cooler than the blackbody. Our 
best fits yielded blackbody contributions of 65-99% to 
the 2-10 keV flux, in stark contrast to the neutron 
star systems where the disk blackbody dominates. The 
blackbody component dominates the 2-10 keV flux most 
when the disk blackbody temperature is lowest. Four 
of our BHCs only exhibited apparent high luminosity 
hard states during 5 ks ACIS observations, and we were 
unable to constrain double thermal models for the re- 
sulting spectra. However, future deep observations may 
strengthen their BHC cases. 

The 31 BHCs that we did fit with the double thermal 
model yielded kT; n 0.6-28<7 outside the NS parameter 
space; the combined \ 2 was 212, for 30 dof (probabil- 
ity ^3xl0 -29 ). Clearly our BHCs are inconsistent with 
being soft NS XBs. 

In Figs. |5]and [7] we plot the 2.0-10 keV luminosity vs. 
temperature for the disk blackbody and blackbody com- 
ponents of the best fit double thermal models for each of 
our BHCs. We overlay the parameter space inhabited by 
the neutron star systems observed by Lin et al. (2007, 
2009, 2012) for spectral states where Comptonization is 
not required in their models, i.e. the atoll soft states 
and Z-source normal and flaring branches. We note that 
Lin et al. quote luminosities from broader energy ranges; 
however, our disk blackbody components make a negligi- 
ble contribution above 10 keV, while the luminosities of 
our blackbody components are typically ~10% higher in 
the 2.0-20 keV band than the 2.0-10 keV band; hence our 
observed 2.0-10 keV luminosity for our BHCs provide a 
fair comparison with the neutron star systems observed 
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Figure 7. Comparison of blackbody 2-10 keV luminosity vs kT 
for the best fit double thermal models to our BHC spectra. Filled 
circles represent field BHCs, and open circles represent GC BHCs. 
The shaded region encompasses the parameter space inhabited by 
hundreds of spectral fits for Aql X-l, 4U 1608-52, XTE J1701-462 
and GX 17+2, which together exhibit every type of NS LMXB 
behavior (Lin et al. 2007, 2009, 2012). We note that the blackbody 
component contributes ~65— 99% of the total 2-10 keV flux for our 
BHCS, but less than 50% for the NS sources studied by Lin et al. 
(2007, 2009, 2012). 

by Lin et al. (2007, 2009, 2012). 

As predicted, our BHCs inhabit a separate region of 
the disk blackbody parameter space to the NS systems 
studied by Lin et al. (2007, 2009, 2012), although some 
sources are consistent with NS values for T m within 3tr. 
We also see a correlation between the 2-10 keV luminos- 
ity and temperature of the disk blackbody component; 
this is due to lower disk temperatures resulting in lesser 
contributions to the 2-10 keV flux. The BHCs are scat- 
tered more widely in the blackbody space, but are also 
systematically offset from the NS systems. 

Remembering that Lin et al. (2009, 2012) failed to 
fit the horizontal branches of XTE Jl 701-462 and GX 
17+2 with their double thermal model, we also mod- 
eled our brightest BHCs (£0.3-10 > 10 38 erg s _1 ) with 
the disk blacbkody + blackbody + Comptonization emis- 
sion mod el that they em ployed for the horizontal branch 
spectra. iLin et all (|2012fl used a cut-off power law to 
model the Comptonization in GX 17+2 over the 2.9- 
60 keV band; they required simultaneous fitting of sev- 
eral spectra to constrain the parameters of the cut-off 
power law, which yielded T — 1.40±0.14, and a cut-off 
energy of 9.9±1.0 keV. We therefore modeled the 0.3-10 
keV XMM-Netwon pn spectra, or 0.3-7.0 keV Chandra 
ACIS spectra of our brightest BHCs with a blackbody, 
disk blackbody and a power law with T fixed to 1.4. 
ILin et all (|2012f ) found that the Compotonized compo- 
nent contributed ~10-50% of the total flux on the hori- 
zontal branch of GX 17+2; we fixed the flux of the power 
law component in our fits to ~30% of the total 2-10 keV 
flux. 

We fitted the spectra of BHCs 1, 3, 16, 20, 21, and 31 in 
this way; the results are summarized in Table |4] We find 
that this three component model results in lower tem- 
peratures for both thermal components than the double 
thermal model fits to the same spectra. I.e., the addition 
of the power law component drives our BHCs away from 
the NS parameter space. Increasing the power law con- 
tribution to 50% of the flux decreased the temperatures 
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Best fit parameters for the blackbody 4 
luminosity, disk blackbody temperature, 
2-10 keV luminosity. Uncertainties are 



Table 3 

disk blackbody model. We provide the absorption , blackbody temperature, 2-10 keV blackbody 
and 2—10 keV disk blackbody luminosity. Finally we give the blackbody contribution to the total 
quoted at the ltr level. The last column states whether the BHC identification is strong (S) or 
plausible (P) 



BHC 


N H 1 10 22 cm" 2 


kT BB / keV 


L BB /10 3 ' 


kT DBB / keV 


iDBB/10 3 ' 


^bb/^tot 


BH 


1 


0.34±0.07 


2 00 +ut> 
z - uu -0.19 


21±2 


0-88_ ; 16 


7.6±0.5 


0.73±0.09 


P 


2 


0.07f 


1 12+ ' 4 
1 - lz -0.11 


3.0±0.5 


50+ ' 3 

u - au -o.io 


0.8±0.4 


0.8±0.2 


P 


3 


0.07f 


i -1O+0.13 
1-1»_0.08 


7.3±0.3 


n 40+0.08 

u - 4y -0.06 


0.8±0.3 


0.90±0.06 


S 


4 


0.19±0.06 


1 50+ 1 ' 
i ' ou -0.19 


3.9±0.8 


0.621°;^ 


0.8±0.4 


0.8±0.2 


P 


5 


0.07f 


I O+25.0 

^-o.s 


2.0±0.7 


n 7 +0.5 
u -'-0.2 


0.6±0.3 


0.8±0.4 


P 


6 


0.07f 


i 34+O.O6 


3.67±0.15 


51 0+ ' 02 
u.oiu_ 013 


0.42±0.05 


0.90±0.05 


s 


7 

8 


0.07f 


1 47+O.O6 
■"-•*' -0.06 


4.59±0.13 


51+003 
u - ol -0.02 


0.54±0.06 


0.89±0.04 


s 
p 


9 


0.07f 


93+ 008 


2.14±0.12 


U.38 o04 


0.14±0.05 


0.94±0.08 


s 


10 


0.07f 


2 4+12.5 


2.3±0.6 


0.86ig; 2 9 


0.7±0.3 


0.8±0.3 


p 


11 


0.07f 


1 11+ - 4 
J - 11 -0.11 


3.0±0.5 


0-50t Q i 


0.8±0.2 


0.80±0.18 


p 


12 


0.07f 


1 24+ ' 7 
1,z -0.19 


1.74±0.17 


59+° 10 
u - oy -0.08 


0.7±0.3 


0.72±0.12 


s 


13 
11 


0.07f 


1 5+ 1 - 4 


1.8±0.4 


fi7+0- 2 
u -°'-0.12 


0.5±0.3 


0.8±0.2 


p 

p 


15 


0.07f 


1 12+ 04 
1 ' 1Z -0.11 


1.6±0.2 


0.56+°;?! 


0.6±0.3 


0.74±0.17 


p 


16 


0.07f 


n n+2.4 
z,z -0.5 


4.5±0.5 


i.03±g:? B 


2.4±0.8 


0.65±0.11 


p 


17 


0.07f 


oq6 +0.13 
u - 9D -0.09 


1.17±0.09 


43+O.O6 
u -*°-o.o5 


0.17±0.07 


0.87±0.10 


s 


18 


0.07f 


1 nc + 0.08 

r-uo_ 05 


4.44±0.19 


41+ ' 06 
u - q:1 -0.04 


0.22±0.08 


0.95±0.06 


s 


19 


0.26±0.04 


1 7+ 1 ' 1 
1- '-Q.3 


3.6±0.5 


0.83+.°;? 9 


1.7±0.5 


0.68±0.13 


p 


20 


0.173±0.010 


1 61+ - 19 
J - Di -0.17 


18.1±1.7 


0.83±S:S 


6.9±1.9 


0.7±0.2 


s 


21 


0.1±0.05 


-i 7+O.8 
11 -0.3 


4.1±0.5 


n 77+0.2 
u -' '-0.13 


1.9±0.7 


0.68±0.12 


p 


22 


0.07f 


i 01+O.O6 

^■^-o.os 


3.71±0.13 


^+ ' 03 
u - oo -0.03 


0.57±0.08 


0.87±0.04 


s 


23 


0.07f 


qqo+0.13 


2.6±0.5 


45+ ' 13 
u -^°-0.06 


0.25±0.13 


0.9±0.2 


s 


21 


0.07f 


i i 9+O.O6 
1 - lz -0.05 


2.12±0.09 


43+003 
u -^°-0.02 


0.19±0.04 


0.92±0.05 


s 


25 


0.07f 


1 5fi+ ' 10 
1 - OD -0.08 


3.22±0.13 


61+ 004 
u - D1 -0.04 


0.39±0.07 


0.89±0.05 


s 


26 


0.07f 


9 q+20.0 
z '°-1.6 


2.1±0.6 


-78±o-.i6 


0.5±0.2 


0.8±0.3 


p 


27 


0.07f 


fl Q4+0- 3 
u - 3 -0.08 


1.56±0.14 


46+ ' 13 
u -*°-0.08 


0.28±0.15 


0.85±0.12 


s 


28 


0.07f 


1 4 +1 '° 


2.3±0.5 


63+ 0,2 
U - D -J_0.12 


0.9±0.4 


0.71±0.2 


p 


29 
30 


0.07f 


1 00+ 007 
1 - uu -0.05 


16.4±0.9 


29+ 006 
u - za -0.04 


0.22±0.11 


0.99±0.08 


s 

p 


31 
32 


0.07f 


1 36+ 005 
1 - OD -0.04 


5.72±0.15 


45+ - 02 
u -*°-0.02 


0.4±0.05 


0.93±0.03 


s 

p 


33 


0.07f 


1 P5+ 03 
1 - lo -0.12 


1.87±0.17 


52+ 012 
u - oz -0.08 


0.38±0.17 


0.83±0.12 


s 


31 


0.07f 


1 5fi+ ' 13 

i - t>D -o.n 


3.7±0.3 


0.63±H 4 4 


0.67±0.12 


0.85±0.08 


s 


35 


0.07f 


1 44+ 1 ' 2 
1 -* 4l -0.19 


2.4±0.4 


0.68±»;i« 


0.8±0.3 


0.75±0.17 


p 



further. Since none of the spectra from our brightest 
BHCs have three component fits that are consistent with 
the fits that successfully described the horizontal branch 
(Lin et al., 2009, 2012), we conclude that none of our 
BHCs are Z-sources. 

4. SUMMARY AND CONCLUSIONS 

We have identified 26 new black hole candidates in the 
central region of M31, using their structure functions or 
luminosities to identify them as X-ray binaries, and their 
high luminosity hard state spectra to classify them as 
BHCs. Of these, 12 are strong candidates, and 14 are 
plausible candidates that may benefit from further obser- 
vations. We were previously limited to identifying BHCs 
in globular clusters, due to the similarities between XB 
and AGN spectra. This brings the total number of BHCs 
within 20' of M31* identified by their high luminosity 
hard states to 35. 

The structure functions of most of our BHCs reveal 
them to be substantially more variable than typical AGN 
(as measured by Vagnetti et al., 2011) over a wide range 



of time scales. Those BHCs with comparable or less 
variability than the AGN have 0.5-10 keV luminosities 
matched b y ^0.6 AGN within the observed region, ac- 
cording to iGeorgakakis etal . (2008). It is therefore un- 
likely that our new BHCs are AGN. 
We have found that our BHC spectra exist in a sep- 

Table 4 

Best fit blackbody and disk blackbody temperatures when our 

brightest BHCs are modeled with a blackbody, a disk blackbody, 

and a Comptonized component represented by a power law with 

T = 1.4 that contributes ~30% of the 2-10 keV flux. This is for 

comparison with the spectral modeling conducted by Lin et al. 

(2012) for the horizontal branch of GX 17+2 



BHC 


N H 1 10 22 cm~ 2 


kT BB / keV 


kT DBB / keV 


1 
3 
16 

20 
21 
31 


0.35 f 
0.07 f 
0.07 f 
0.197±0.011 
0.11±0.07 
0.07 f 


1 oo+i-«E> 

1 - yy -0.4 
1 07+ ' 19 

l-44±o.3 

1 29+ ' 2 
i - ZM -0.13 

l-6±o.4 

-, 97 +0.09 
-"-■ z '-0.08 


i+l.l 
*--$.?, 

0.42+°f 

o.8i°;'t 

0.73±°;il 

7 +0 - 4 
u -'-0.2 

n 43+003 

u - 4J -0.03 
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arate parameter space to Galactic neutron star systems 
when we compare our best fits for absorbed blackbody + 
disk blackbody emission models with the systems stud- 
ied by Lin et al. (2007, 2009, 2012). This is expected 
because this double thermal model fails to fit hard state 
spectra; the disk blackbody component is forced low tem- 
peratures in order to provide the low energy flux. Indeed, 
the probability that our BHCs all lie in the NS parameter 
space when fitted with the double-thermal model is just 
~3x 10~ 29 . Seven of our BHCs (20% of our total sample) 
were found within 100" of the M31 nucleus, lending sup- 
port to the hypothesis that the M31 bulge is sufficiently 
dense to form a significant number of XBs dynamically, 
as seen in globular clusters; since the stellar velocities in 
the M31 bulge are considerably higher than in GCs, sur- 
viving XBs are expected to have short periods and are 
likely to contain black hole accretors (jVoss fc Gilfanovl 



likely 
[2001 . 
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BHC7 N H = 1 .6+/-0.2 (34/85) T = 1 .69+/-0.05 (44/85) BHC8 N„ = 1 .5+/- 1 . 1 (2/4) r = 1 .4+/-0.2 (2/4) 
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BHC 1 3 N = 1 .2+/-0.3 (5/8) r = 1 .74+/-0.09 (5/8) 



BHC14 N H = 4+/-2 (1.5/5) r = 2.6+/-0.5 (2/5) 
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BHC19 N„ = 5.0+/-0.6 (16/18) r = 


1.85+/-0.09 (11/18) 
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BHC21 N„ = 2.1+/-0.2(27/94)r= 1.71+/-0.04 (40/94) 





Log(x/day) 
BHC22 N H = 1.5+/-0.2 (22/66) r = 1 .67+/-0.05 (42/66) 
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BHC24N H = 1.8+/-0.4(3/15)r= 1.83+/-0.18 (3/15) 




Log(T/day) 



Log(x/day) 



Figure 5. continued 
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BHC25N =0.9+/-0.3(6/16)r= 1.32+/-0.06(11/16) 




BHC26N H = 1.6+/-0.4(9/17)r= 1.98+/-0.10 (7/17) 



BHC29N H = 0.8+/-0.7 (3/1) r= 1.46+/-0.13 (18/1) 
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BHC28N H = 1.6+/-0.4(3/18)r= 1.78+/-0.09 (13/18) 
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BHC30 N H = 2.9+/-1.5 (1.4/4) r = 1.9+/-0.3 (3/4) 
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Figure 5. continued 
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Barnard et al. 

BHC31 N H = 2.34/-0.17 (24/84) r = 1.50+/-0.03 (35/54) BHC32 N„ = 1.2+/-1.2 (0.1/2) T = 1.6+/-0.3 (0.1/2) 




Log(T/day) 
BHC35 N H = 2.1+/-0.9 (3/7) T = 1.81+/-0.17 (7/7) 
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BHC34 N H =1.3+/-0.4 (3/16) r = 1.51+/-0.09 (8/16) 
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Figure 5. continued 



